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To understand the non-exponential relaxation associated with solvation dynamics in the ionic 
liquid l-ethyl-3-methylimidazolium hexafluorophosphate, we study power spectra of the fluctuating 
Franck-Condon energy gap of a diatomic probe solute via molecular dynamics simulations. Results 
show l//dependence in a wide frequency range over 2 to 3 decades, indicating distributed relaxation 
times. We analyze the memory function and solvation time in the framework of the generalized 
Langevin equation using a simple model description for the power spectrum. It is found that the 
crossover frequency toward the white noise plateau is directly related to the time scale for the 
memory function and thus the solvation time. Specifically, the low crossover frequency observed in 
the ionic liquid leads to a slowly-decaying tail in its memory function and long solvation time. By 
contrast, acetonitrile characterized by a high crossover frequency and (near) absence of 1// behavior 
in its power spectra shows fast relaxation of the memory function and single-exponential decay of 
solvation dynamics in the long-time regime. 



I. INTRODUCTION 

Solvation dynamics in room-temperature ionic liquids 
(RTILs) have received intensive theoreticalii^i^i^i^i^i^'^ 
and experimental attention recently .^^'^^i^^'^'^i-'^'^'^^i^^ The 
collective influence of the solvent in the presence of an 
optically-active probe solute is usually monitored via var- 
ious dynamic electronic spectroscopies and described in 
terms of time-dependent fluctuations and relaxation of 
the Franck-Condon (FC) transition energy of the solute 
in solution. Often observed in RTILs is biphasic re- 
laxation comprised of ultrafast sub-picosccond dynamics 
and ensuing non-exponential decay. Short-time solvation 
dynamics, arising mainly from small-amplitude inertial 
translational motions of solvent ions, make a substantial 
contribution to overall solvent relaxation despite their 
high viscosity,-^ and thus can play an important role 
in reaction dynamics^^ in RTILs. The subsequent re- 
laxation, attributed to diffusive dynamics, involves coop- 
erative movement of ions and accompanying structural 
relaxation. 

A spectral analysis of the FC energy gap, which is 
equivalent to analyzing its time correlation function, 
is useful for investigating dynamics over various time 
scales. Of particular interest are long-time fluctuations 
that often result in non-exponential relaxation. For in- 
stance, 1// behavior for long-time fluctuations in water^^ 
and liquid silica— has been reported. Defect fluctu- 
ations in a disordered two-dimensional liquid also ex- 
hibit 1// spectra, suggesting that system dynamics are 
hctcrogeneouSi^ Slow relaxation in these systems indi- 
cates that the structural memory persists for a long time 
in spite of the fast inertial motion of solvent molecules. 

In the generalized Langevin equation (GLE) descrip- 
tion of solvation dynamicsj^ii^ the time-dependent fric- 
tion plays the role of a memory function. The tem- 



poral behavior of the memory function for a variety of 
liquid systems has been analyzed via numerical trans- 
forms of appropriate time correlation functions^ i^^'^'^i^^ 
and a self-consistency method— The memory function 
couched in terms of radial distribution functions has also 
been studied numerically.^^ Memory functions thus ana- 
lyzed have two general characteristics in common, a rapid 
decay within the first few hundred femtoseconds of relax- 
ation and a long tail thereafter. The former among these 
features originates mainly from short-time collisions via 
fast inertial motions of molecules in liquids. The residual 
memory effect which can last for a prolonged period gov- 
erns the long-time behavior of relaxation dynamics. In 
solvation dynamics, the overall friction given by the in- 
tegration of the memory function including its long-time 
tail is directly related to the solvation time. 

In this study, we investigate the power spectra of 
the FC energy gap in the room-temperature ionic liq- 
uid, l-ethyl-3-methylimidazolium hexafluorophosphate 
(EMT'^FF^), employing a diatomic probe solute via 
molecular dynamics (MD) computer simulations. Among 
others, 1// dependence and crossover to the white noise 
are observed in the low frequency region of the spectra. 
The latter represents the onset of normal diffusion, which 
marks the escape from the subdiffusive regime in the in- 
termediate time scale corresponding to the 1// region— 
In the GLE description of solvation dynamics, the mem- 
ory function can be expressed conveniently in terms of the 
power spectrum of the FC energy gap. Using a simple 
model description for the power spectrum, we compute 
the memory function and confirm that the crossover fre- 
quency determines the time scale of the memory function 
decay and is inversely proportional to the solvation time. 

This paper is organized as follows: In Sec. [Til we 
briefly review solvation dynamics, its GLE description 
and power spectrum analysis. In Sec. Illll simulation 
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methods are described. MD results for the time cor- 
relation functions and power spectra are presented in 
Sec. IIV| while the memory functions are analyzed by us- 
ing a model description in Sec. [V] Concluding remarks 
are offered in Sec. IVll 



we obtain 



II. SOLVATION DYNAMICS 

In this section we review briefly the time correlation 
function of the FC energy gap and its power spectrum, 
and the GLE approach to equilibrium solvation dynam- 
ics. We assume that the probe solute is characterized 
by an active electronic state a and a reference electronic 
state r as in previous MD studies.— i^i^ For a given solvent 
configuration, the FC energy associated with the a ^ r 
transition of the solute is given by 



Er — Ea 



(1) 



where Ea^r denotes the total energy of the solute-solvent 
system with the solute in state a and r, respectively. 

We consider the equilibrium solvation dynamics, char- 
acterized by the normalized time correlation function 

^'^/'■^'^ = mE^^rY) ' 

where (• • • ) denotes the equilibrium ensemble average 
in the presence of the a-state solute and 5lS.Ea^r = 
AEa^r — i^Ea^r) depending on time t. Henceforth the 
subscripts representing the solute electronic states are 
suppressed for brevity. 

The time correlation function is conveniently de- 
scribed by the GLE,^-'^'^^ derived via the Mori-Zwanzig 
projection?^'^^ onto a set of dynamical variables 
{SAE,SAE}: 

SAE{t) = -ojI SAE{t) - [ dt' C{t - t') SAE{t') + R{t), 

Jo 

(3) 

where the solvent frequency lUs — 



((5Aij)2)(((5A£^)2)-i characterizes inertial dy- 
namics of A£^j^- and the time-dependent friction (i.e., 
the memory function) and the random force R{t) 
(scaled by the inertia associated with SAE dynamics) 
are related via the fluctuation-dissipation theorem 



(i?(0)i?(i)) 



(4) 



{{SAEr) 

The equilibrium time-correlation function C{t) satisfies 



C{t) = -ujlC{t) - / dt'C{t - t')C{t') . 
Jo 

Taking the Laplace transform 



C(z) 



exp(-zt) C{t)dt , 



(5) 



(6) 



C{z) 



(7) 



which shows that the solvation time Tgoiv = /q°° C{t)dt = 
(5(0) is governed by C(0) and cj^. 

The normalized power spectrum of the FC energy gap 
SAEa^r is defined to be 



Sa/rif) 



lim — 



J-T/2 ^{{SAEa^rY) 



(8) 

which corresponds to the Fourier transform of the (nor- 
malized) correlation function C(t), according to the well- 
known Wicncr-Khintchine theorem.— Namely, the corre- 
lation function, which is real, relates to the power spec- 
trum via 



C{t) = 2tt / dfS{f)cos2Trft 



(9) 



Because of the discrete nature of sampling in simulations, 
we use the discretized form of the power spectrum!^ 



Sifk) 



At 

'n' 



SAE.„ 



^2-iTink/N 



(10) 



where N is the number of samples, fk (= k/NAt) is the 
fcth frequency (fc = 0, ■ • ■ , A''— 1) and At is the sampling 
interval. Eq. ^ shows that the power spectrum S{f) 
carries all the information relevant to the time correlation 
function C (t) and elucidates its behavior in the frequency 
domain. Therefore S{f) can reveal insight into dynamic 
character and time scales of the system. For example, in 
the case of Debye relaxation given by single-exponential 
decay, S{f) is a Lorentzian that decreases as 1//^ at high 
frequencies. If the system is characterized by multiple 
relaxation times, say, f^^ and /j"^ (with /i <^ /2), its 
power spectrum displays a plateau in the low frequency 
region (/ <^ fx) and 1//^ decay in the high frequency 
limit f ^ f2- In the intermediate region /i ^ / ^ /2, 
S{f) varies approximately as l//<^ 



III. SIMULATION METHODS 

The simulation cell is comprised of a rigid diatomic 
solute immersed in EMI^PFg", consisting of 112 pairs of 
rigid cations and anions. We consider two different solute 
charge distributions: a neutral pair (NP) without charges 
and an ion pair (IP) with unit charge separation. When 
solvent dynamics occurs in the presence of NP, i.e. NP is 
the active electronic state, IP is regarded as the reference 
electronic state, and vice versa. The Lennard-Jones (LJ) 
parameters employed for each constituent atom of the 
solute are a = iA and e/ks = 100 K (with the Boltz- 
mann constant ks)- Its bond length remains fixed at 
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3.5 A for all cases considered here. As for the solvent in- 
teraction potential, we employ the same parametrization 
as in Ref . l33l 

MD simulations were conducted in the canonical en- 
semble at temperature T = 400 K through the use of 
the DL_POLY programi^ For each solute charge distri- 
bution, we simulated the combined solute-solvent sys- 
tem for 70 ns after 6 ns equilibration. To compute the 
power spectrum of the energy gap with reduced noise, 
we divided the trajectory into seven 10 ns segments and 
calculated <S'(/) in Eq. pH)) by averaging over the seven 
segments. The sampling time interval At in our simula- 
tions was 10 fs, which sets the maximum frequency of our 
power spectrum analysis at 50ps~^, while the minimum 
is 10-4 ps-i. 

For comparison, we also performed simulations in apro- 
tic acetonitrile at T = 300 K. In the simulation cell, a sin- 
gle NP (or IP) solute is immersed in 512 rigid molecules 
of acetonitrile. The LJ parameters and partial charges 
for acetonitrile were taken from Ref. |3^ The trajectory 
was 20 ns long and the FC energy gap was saved at every 
2fs. 



IV. TIME CORRELATION FUNCTION AND 
POWER SPECTRUM 

In this section, we present MD results for equilibrium 
solvation dynamics in EMI^PFf^ and acetonitrile. We 
begin with the time correlation functions of the FC en- 
ergy gap fluctuations in the presence of the NP and IP 
solutes in Fig. [T] 

The MD results for Ca/rit) in EMI+PFg" exhibit 
biphasic relaxation, i.e., ultrafast inertial relaxation fol- 
lowed by an extremely slow decay [Fig.[TJa) and (b)], con- 
sonant with prior simulation studies ji'^ii^ii^i^'^ Long-time 
behaviors of Ca/r{t) are well described by the stretched 
exponential function exp[— (t/ro)^]:— "^"-^ A good agree- 
ment with the MD results was obtained with fitting pa- 
rameters /3 = 0.30 and tq = 1.75 ps for Ips < t < 50ps 
in the case of NP. For IP, the stretched exponential fit 
with f3 = 0.18 and tq = 16.84 ps apphes to a time range, 
Ips i$ i i$ 1000 ps, which is about 20 times wider than 
the NP case. The result that the tq value with IP is 
larger than that with NP by one order of magnitude is 
a direct consequence of slow solvent relaxation dynamics 
in the presence of the former solute, compared with the 
latter. For instance, the Ca/r{t) value reduces to below 
0.1 after t > Ins in the presence of IP, whereas it takes 
much shorter ^ 34 ps with NP. The corresponding solva- 
tion times are Tgoiv = 450 and 1 1 ps for the IP and NP 
solutes, respectively. 

For comparison, we consider Ca/rit) in acetonitrile in 
Fig. [Ijc) and (d). We notice that solvation dynamics 
in acetonitrile are much faster than those in EMI^PFg", 
congruent with previous studies. The MD results for 
Tsoiv for NP and IP are 0.17 and 0.29 ps, which are 
smaller than the EMI~''PF^ values by 2-3 orders of mag- 



nitude. The biexponential functions cexp{—t/ti) + (1 — 
c)exp(-</<2) with c = 0.97, h = 0.14 ps, ^2 = 1.086 ps 
for NP and c = 0.85, h = 0.09 ps, ^2 = 1.31 ps for 
IP provide excellent fits for Ca/r{t) in acetonitrile for 
t ^ Ips- As analyzed below, Sa/r{f) associated with 
Ca/r{t) in acetonitrile does not show pronounced 1// 
behavior. Thus while a stretched exponential function 
also appears to yield a reasonable fit (result not shown 
here), we will take the view that Ca/rit) in acetonitrile 
is biexponential with two relaxation times. In the long- 
time limit of this description, Ca/rit) becomes a single 
exponential decay characterized by the longer of the two 
relaxation times. 

To gain additional insight into characteristics of equi- 
librium solvation dynamics, we consider the power spec- 
trum Sa/rif) of the fluctuating FC energy gap AEa^rit)- 
Figure [2] presents Sa/rif) in EMI^PFg" and in acetoni- 
trile, obtained from MD simulations with the aid of 
Eq. (Uni). We observe in Fig. ^a) that Sa/rif) of NP 
in EMI^PF^ is characterized by at least four differ- 
ent regimes >i^ii2ii2^ In the low frequency region below 
0.001 ps~^, the spectrum is of a white-noise type, which 
reveals the absence of correlation between two events sep- 
arated by longer than ~ 1 ns. For later use, we denote 
as fc i "crossover frequency" ) the frequency below which 
Sa/rif) flattens out. For the NP solute under consid- 
eration here, fc ~ 10~^ps~^. As / increases above fc, 
the power spectrum exhibits the 1// behavior, which is 
often interpreted as the presence of many different re- 
laxation time scalesi^ If the frequency further increases 
beyond ^ lps~^, Sa/rif) begins to drop sharply. This 
is attributed to the rapid initial decay of Ca/rit) at 
short times [cf. Fig.[lja)]. Finally, in the high frequency 
region / > 10ps~^, the power spectrum decreases as 
^ 1//^. For perspective, the reader is reminded that 
/ = 10~'*ps~^ and 50ps~^ are the lower and upper lim- 
its of our analysis of the MD results in the frequency 
domain. 

The power spectrum for IP in EMI+PFg in Fig.H^b), 
though similar to that for NP, shows a couple of in- 
teresting differences. First, we were not able to ob- 
serve a white-noise spectrum in the low frequency re- 
gion of Sa/rif) in the presence of IP. The most hkely 
reason is that the MD trajectory is not long enough to 
probe the onset of the complete loss of correlations be- 
cause Ca/rit) in the IP case decays much more slowly 
than that in the NP case. We ascribe this difference in 
Ca/rit) to electrostriction. To be specific, it tends to 
enhance the solvation structure around IP compared to 
NP and in turn makes rugged the landscape of the po- 
tential energy surface, upon which solvent ions diffuscii 
This leads to slower diffusion of solvent ions and there- 
fore slower relaxation of AE fluctuations in the pres- 
ence of IP than NP. We thus expect that the white-noise 
plateau for the former solute will appear at frequencies 
lower than 10~^ps~^. This means that fc for IP would 
be lower than fc for NP, 10~^ps~^, by more than one 
decade. Second, Sa/rif) with IP shows a flat region for 
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O.lps-i < / < Ips-i. While Sa/rif) with NP also has 
a hint of a plateau around / = lps~^, it is much more 
prominent in the case of IP. The plateau is a part of the 
Gaussian power spectrum associated with the ultrafast 
initial relaxation of Ca/r{t). One quick and easy way 
to see this is that infinitely fast delta function relaxation 
yields a plateau in the frequency domain. This plateau is 
more noticeable in IP because the difference in time scales 
between its inertial and diffusive dynamics is larger than 
that of NP. 



The power spectra of the stretched exponential fits in 
EMI+PFg are compared with the MD results in Fig.H^a) 
and (b). For both IP and NP, the stretched exponen- 
tial functions well describe the 1// characteristics — both 
the range and exponent of the power law behavior — 
of the simulation results. One prevalent notion is that 
stretched exponential behavior arises from the superpo- 
sition of different single exponential decays, weighted by 
a broad distribution of relaxation timesi^ While the as- 
sumption of single exponentials per se may be too re- 
strictive, there is considerable evidence that RTIL dy- 
namics are characterized by a distribution of different 
time scales i^i^^i^i^i^ As pointed out in Sec. HH the 
power spectrum of the system, involving a number of 
relaxation processes of different time scales, in general 
exhibits 1// dependencci^ Thus the presence of the 1/f- 
type domain in the power spectrum of AE is another 
manifestation of non-exponential RTIL relaxation. This 
also suggests the similarity between solvation dynamics 
in the RTILs and the dynamic heterogeneity observed in 
glassy liquids 



Before we turn to memory functions, we briefly con- 
sider Sa/rif) for acetonitrile in Figure [D^c) and (d). In 
the case of NP, Sa/rif) displays a white- noise plateau 
for / ^ 2ps^^ but it does not show a 1// power-law 
behavior. For IP, the white-noise spectrum obtains for 
/ < 0.08 ps~^. As / increases beyond the white-noise re- 
gion, Sa/rif) briefly shows 1//" behavior with a ~ 0.53 
before it begins a rapid decrease. As in EMI^PFg", the 
crossover frequency fc for IP is smaller than that for 
NP, again due to electrostriction. Regardless of the so- 
lute charge distributions, the onset of the white-noise re- 
gion is at much higher frequencies in acetonitrile than 
in EMI^PFg because the long-time solvent relaxation in 
the former is much faster than that in the latter. Fur- 
thermore, the 1// spectrum is nearly absent in Sa/rif) 
of acetonitrile in contrast to the EMI^PFg' case. As 
mentioned above, this is why we favor a biexponential 
description for Ca/rit) hi acetonitrile over a stretched 
exponential description. We notice that except for the 
rapid decrease in the neighborhood of / = 10ps~^ arising 
from ultrafast inertial relaxation of Ca/rit): the biexpo- 
nential functions indeed provide an excellent framework 
to describe Sa/rif) acetonitrile. 



MEMORY FUNCTION ANALYSIS 



Here we consider a simple model description for 
and analyze its memory function ^(t) to gain insight into 
long-time solvation dynamics. As observed above, power 
spectra in EMI~''PF^ are characterized by at least four 
different regimes; viz., Sif) generally show white-noise, 
1//, sharp fall-off and 1//^ behaviors in turn as the 
frequency increases. Since we are mainly interested in 
long-time dynamics, we ignore the frequency region (i.e., 
Ips""^ ^ / ^ 10ps~^) of rapid Sif) decay, which arises 
mainly from ultrafast inertial relaxation of C it) . We thus 
consider the power spectrum of the following form: 



Sico) 



Auj~^ for < w < ojc , 

Auj~^ for LUc < < (^0 

AuJqLU^'^ for LUq < UJ < UJm 

for > LOm , 



(11) 



where Si— to) = Siuj), lo = 27r/ and A is the normaliza- 
tion constant given by 2 A = [2 — (jjo/f^m + lii(wo/wc)] 
By differentiating Eq. ([9]) twice with respect to t and 
setting i = 0, we derive a sum rule^^ 



duj uP'SiuS) 



2 - ujo/uom + ln(wo/^^c) ' 



(12) 



which clearly shows the necessity of a high frequency cut- 
off at LJrn in Eq- (fTT|) . 

The parameters employed to model >S'(/) are compiled 
in Table HI SI and S2 there are the model descriptions 
for solvation of NP and IP in EMI'^'PF^, respectively, 
while the corresponding cases in acetonitrile are modeled 
by S3 and S4. Thus SI to S4 correspond to the cases in 
Fig. [2ja) to (d), respectively. The solvent frequency los 
was obtained from the simulations via Eq. ([3]) and the 
cutoff frequency ujm was estimated from Eq. p2)) . The 
values of lOc and ojq were determined with the aid of the 
MD resuhs in Fig. El In the presence of IP in EMI+PF^, 
we were not able to obtain lOc from the simulation because 
the plateau behavior in the low frequency region of Sif) 
was not accessible as mentioned above [Fig.jSfb)]. In the 
absence of any additional information on the crossover 
frequency, we assumed uj^ — 21: x 10^* ps^^, the lowest 
frequency that is allowed in our analysis of MD results, 
for S2. The corresponding crossover frequency for SI, 
viz., NP in EMI+PFg", is lo^k 2ttx IQ'^ ps'^ For S3, we 
employed lOc = uja because acetonitrile does not exhibit 
a 1// spectrum in the presence of NP. The narrow 1// 
region present in S4 was incorporated into the model cal- 
culations by choosing ujc — 0.167rps~^ and cjq = 47rps~^. 

Within the model description of Eq. (|lip . we can de- 
termine exactly. Specifically, we make Laplace trans- 
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form of Eq. (O with Eq. ([TT]) to find 



Ciz) 



z 



= 2A 



US'' + 

z 



duj- , „ „, 
UJc[U)^ + z^) 



duj- 



LUqZ 



= 2A 



UJ^{UJ^ + z^) 



— tan"^ — H In- 



UJc 



2z 



duj 



^0 



(jj{u!'^ + Z^) 



^0 



-In- 



Lui + Z^ 



Wo 
"T2 



ftan^i 
z^ \ 

We rewrite Eq. ([7]) as 



Z Vti^o 



tan 



1 - zC{z) 



(13) 



(14) 



substitute Eq. (|13p into Eq. and make inverse 

Laplace transform^^ numerically to obtain C(t). 

Figure [Slja) displays the results for the memory func- 
tion C(i) thus obtained. For comparison, we employ 
the method used in Ref. to determine directly 
from C{£) and associated time correlation functions of 
a nonconservative force'^^ and present the results in Fig- 
ure [3]Jb). We notice in Figure El^a) that C(i) falls off 
very quickly at the very early stage of relaxation, regard- 
less of the crossover frequency Wc- Subsequent decay of 
the residual memory varies strongly with Wc- Generally, 
the long-time relaxation of C(t) becomes slower with de- 
creasing u)c^ Comparison of the results in Figure [3]^ a) 
and (b) shows that while there are differences, C(t) ob- 
tained from model S{lS) in Eq. (jlip correctly captures im- 
portant features of memory both at the qualitative and 
semi-quantitative level. These include rapid initial de- 
cay of a large amplitude and the existence of a long-time 
tail. Nonetheless, the model calculations yield spurious 
oscillations at short times (see below) and overestimation 
of the long-time memory effect. For example, the model 
predictions for the magnitude of the residual memory, say 
at t ~ Ips, are considerably larger than the MD results. 
This overestimation is attributed mainly to the neglect 
of the rapid decay of S'(/) in our model description at 
high frequencies. 

Here we briefly examine the oscillatory behavior of 
C(i) observed in the first few hundred femtoseconds in 
Fig. [3)Ja). It is of interest to note that the period of 
these oscillations is close to 2Ti/u)m- This suggests that 
oscillations are closely linked to the presence of a cutoff 
at frequency in the model power spectrum employed 
in this study. To check this, we considered a Gaussian 
power spectrum S{f) = {A/ujq) exp[— (w — cjo)^/2Wg] for 
Lo > loq that does not require a high-frequency cutoff. 
The resulting friction is exhibited in Fig. 21 The dis- 
appearance of the rapid initial oscillations confirms that 



indeed the cutoff is mainly responsible for rapid oscilla- 
tions in C,{t) in Fig. [3{a). We note that slow oscillations 
in Fig. U arise from the discontinuity in the derivative of 
the model Gaussian spectrum aX uj = ujq. 

Finally, we consider solvation time Tso1v(= C'(O))- It is 
related to the total memory C(0) via Eq. (O, so that 



C(o) 



'^solv 



TT 
2uJc 



WQ 



In 



(15) 



where we have used Eq. (|13p in passage to the final ex- 
pression. The results for Tgoiv obtained from Eq. (jlSp are 
presented in Table |T1 We notice that TsoIv is inversely 
proportional to the crossover frequency uJc- Thus, all 
other things being equal, TsoIv increases as uJc decreases. 
This is directly related to the observation made above 
that long-time relaxation of C^{t) becomes slower with de- 
creasing ojc- Also interesting is that TsoIv varies with the 
frequency range I^c associated with the 1// behavior 
in S{f). For instance, for given lUc, the growing range 
of the l/f region, i.e., increasing loq, tends to reduce 
the solvation time. Since the lower limit of distributed 
time scales for l/f is ujq^ [cf. Eq. pTjl ]. the shortest 
time scale for l/f becomes faster as ojq increases. In 
other words, faster processes become available increas- 
ingly more to the system, while the availability of slower 
processes remains unchanged. This yields the reduction 
in the solvation time because contributions to solvent re- 
laxation from faster processes become progressively more 
important than slower processes. 



VI. CONCLUSIONS 

We have studied the time correlation functions and 
power spectra of the FC energy gap and related memory 
functions associated with solvation dynamics of model 
diatomic solutes in EMI+PFg" and CH3CN. It was 
found that the power spectra of both the NP and IP 
solutes in EMI"'"PFg display l/f dependence over a 
range of intermediate frequencies, which indicates non- 
exponential relaxation dynamics in RTILs. In the case 
of NP, we observed white-noise behavior at low frequen- 
cies. The power spectrum of NP in acetonitrile does not 
yield the l/f dependence and its crossover frequency is 
much higher than the corresponding value in EMI^PF^. 
Though not pronounced, the IP solute in acetonitrile 
shows l/f dependence but in a frequency range much 
narrower than that in EMI"'"PF^. 

Using a simple but analytic model for power spectra, 
we have determined the memory function in the GLE de- 
scription of solvation dynamics and compared with the 
MD results. With proper account of white-noise, l/f 
and 1//^ behaviors of power spectra, we have found that 
the model description reproduces the MD results of fric- 
tion reasonably well. We have also found that the time 
scale of memory effects and the l/f regime are closely 
related. We have obtained the solvation time in terms of 
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the frequency parameters of the model power spectrum 
description. Among others, the solvation time was found 
to be inversely proportional to the crossover frequency. 
Together with our MD results of Sa/r{.f)i this indicates 
that the solvation time for the IP solute in EMI^PF^ is 
far longer than that for NP, at least by one order of mag- 
nitude. Electrostriction which exerts a strong influence 
on the landscape of the potential energy surface relevant 
to solvent ion diffusion is mainly responsible for the varia- 
tion of solvation time with the solute charge distribution. 
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In the numerical inversion of Eq. (14) to the time domain, 
the linear term in z should be cancelled out by asymptotic 



behavior of lij1C{z) l\i — C(z)] for large z. Unlike the an- 
alytic expression in Eq. (13), however, rmnierical results 
of C{z) do not satisfy this condition. We have thus de- 
termined C,[z) by combining small and large z behaviors 
obtained from C(t) and from the time correlation function 
of the nonconscrvativo force, respectively. 
We have performed additional calculations by varying lOc 
with Wo and fixed and confirmed that C(t) decays more 
slowly with decreasing Wc- 
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TABLE I: Frequency parameters for the power spectrum in Eq. (|11|) and the solvation time TsoIv given by Eq. (|15|) . The 
frequency and time are measured in the units of ps"^ and ps, respectively. SI through S4 model the power spectra obtained 
from MD simulations of solute-solvent systems shown in the last column. 




FIG. 1: Time correlation function C„/^(t) of (5AK^r(t) in (a) EMI+PFg for NP/IP, (b) EMI+PFg for IP/NP, (c) CH3CN for 
NP/IP, and (d) CH3CN for IP/NP active/reference states of the solute. Insets in (a) and (b) display semi-log plots for Ca/rit) 
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FIG. 2: Power spectrum Sa/r{f) of 5AEa-tr{t). The results obtained via Eq. (I10|l are plotted in "+" symbols, while the Fourier 
transforms of Ca/r(i) are given in a solid line. They show an excellent agreement as they should. Also displayed are Fourier 
transforms of the stretched exponential and biexponential fits, (a), (b), (c), and (d) refer to the same cases as in Fig. [1] 
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FIG. 3: Memory function (^{t), (a) obtained via the inverse transform of Eq. (|14p . together with Eq. (|13|l and the parameters 
for SI to S4 in Table U (b) evaluated directly from the correlation function for the NP and IP solutes in EMI"'"PF^ and in 
acetonitrile. 




FIG. 4: Memory function (^{t) in case that the power spectrum at high frequencies (u > ljo) follows a Gaussian function: 

S{ijj) — (A/loo) exp[—{uj — u!o)'^/2iOg] with ujq = 27rps~^ and utg set to be O.Swo- The solvent frequency is given by lOb = 5.7ps~'^ 

-3 , 



for cjc = 27r X 10 ps ^ according to the sum rule. 



